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Решение задач расчета прочности 
конструкций на МИМО-компьютере 


Рассматривается решение задач расчета прочности конструкций на компьютерах с параллельной 
организацией вычислений, в том числе алгоритмы формирования, исследования и решения систем 
разрешающих уравнений, получаемых при дискретизации методом конечных элементов названных 
задач. Приводятся некоторые результаты численных экспериментов на МГМО-компьютере. 


Введение 


При проектировании объектов различного назначения часто возникает необходи- 
мость в проведении расчетов прочности конструкции таких объектов или отдельных 
ее элементов. В первую очередь задачи расчета прочности конструкций возникают в 
различных отраслях машиностроения (судостроении, авиастроении, ракетостроении, 
моторостроении и т.п.), а также в промышленном и гражданском строительстве 
(расчет отдельных конструктивных элементов или сооружений в целом). 

Для проведения прочностных расчетов создано множество программных средств, 
как универсального назначения, так и ориентированных на определенную отрасль 
(МАЗТКАМ, МАКС, АМЗУЗ и другие, см. также [1], [2]). Например, для решения 
задач расчета прочности строительных конструкций разработан программный 
комплекс (ПК) ЛИРА [2], [3]. Подстроив его интерфейс для соответствующей 
предметной области, ПК ЛИРА можно использовать для прочностных расчетов в 
других отраслях. 

Возрастающие требования к качеству проектных решений, применение новых 
конструктивных материалов вызывает необходимость в решении качественно новых 
задач, а также в выполнении расчетов уникальных конструкций. Растет потребность 
в новых методах и подходах, связанных с построением и исследованием корректных 
компьютерных моделей, которые адекватно отображают реальную работу конструкций. 
Возрастают требования к достоверности получаемых компьютерных результатов. 
Эти факторы приводят к тому, что существенно возрастает объем обрабатываемой 
информации, что, в свою очередь, вызывает рост требований к вычислительным 
ресурсам. 

Разрешить эту проблему позволят использование в прочностных расчетах МИМО- 
компьютеров и интеллектуализации процессов исследования и решения задач. Так, на 
базе интеллектуальных программ библиотеки шра® [4], [5], препроцессора и постпро- 
цессора ПК ЛИРА разработан программный комплекс по исследованию и решению 
некоторых задач расчета прочности строительных конструкций на интеллектуаль- 
ных параллельных компьютерах Инпарком [6]. 
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Математические постановки задач 


Математически задачи расчета прочности конструкций с использованием принципа 
возможных перемещений могут быть поставлены в виде следующих вариационных 
задач [7]: необходимо найти вектор-функцию ие, которая для любой вектор- 
функции уЕЦу (любого возможного перемещения) удовлетворяет одному из интег- 
ральных тождеств: 

— для статической задачи 


а(и,у) = Ку); (1) 
— для динамической задачи 
а(и,у) + Б(и"у) + с(и’») = КУ, ию=и®, и’) =ио; (2) 
— для задачи о собственных колебаниях 
а(и,») = АБ(и,»), (3) 


где Ци — бесконечномерное функциональное пространство возможных перемещений, 
симметричные билинейные функционалы а(и,у), Б(и”у), с(и’у) пропорциональны 
соответственно потенциальной, кинетической энергиям деформации и энергии тормо- 
жения, а линейный функционал (ру) пропорционален работе приложенных (внешних) 
усилий при нагружении. Через и’ обозначена первая производная вектор-функции 
и(1,)) по времени, а через и ”- вторая. 

Здесь рассмотрены только линейные задачи, поскольку предполагается, что 
решение нелинейной задачи может быть сведено к решению последовательности ли- 
нейных задач. 


Дискретизация задач методом конечных элементов 


Теоретической основой большинства программных средств для расчета проч- 
ности конструкций является метод конечных элементов (МКЭ) [8], реализованный в 
форме перемещений. Выбор именно этой формы объясняется простотой ее алгорит- 
мизации и физической интерпретации, наличием единых методов построения матриц 
жесткости и векторов нагрузок для различных типов конечных элементов, возмож- 
ностью учета произвольных граничных условий и сложной геометрии рассчитываемой 
конструкции. 

Дискретизация вариационной задачи методом конечных элементов заключается 
в замене бесконечномерного пространства допустимых функций (4 его конечномерным 
подпространством Со’ кусочно-полиномиальных функций (пространством метода 
конечных элементов). Для построения такого пространства исходная область разби- 
вается на подобласти-элементы. На каждом элементе функции из пространства МКЭ 
являются полиномами заданной степени, коэффициенты которых определяются через 
значение функций (а в ряде случаев и их производных) в узлах элемента (точки в 
вершинах, на сторонах, гранях или внутри элемента). Таким образом, функции тако- 
го пространства МКЭ однозначно определяются конечным числом своих значений 
(и значений своих производных) в узлах элементов, которые образуют в исходной 
области задачи, включая границу, конечно-элементную сетку. Размерность прост- 
ранства МКЭ | равна общему количеству этих неизвестных значений в узлах сетки 
(эти значения называют еще узловыми параметрами или степенями свободы). 
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Узловым параметрам соответствуют базисные (координатные) вектор-функции 
фь отличные от нуля только на элементах, которые содержат соответствующий дан- 
ному узловому параметру узел. Для узловых параметров и базисных функций 
справедливы следующие соотношения 


ГАФ) = 6, (4) 
где дб, — символ Кронекера, а результатом операции Ё(ф;) является значение компо- 
ненты вектор-функции ф; для степени свободы (узлового параметра) Г,. Вектор- 
функции из подпространства (' могут быть представлены в виде удовлетворяющей 
главным (кинетическим) краевым условиям линейной комбинации базисных вектор- 
функций 


и" (ч) = Уха (ч), (5) 


где ф; ((=1,2, ..., п) - упомянутый выше кусочно-полиномиальный базис Оо". 

Тогда билинейные и линейный функционалы из (1) - (3) в конечномерном про- 
странстве Оо’ можно записать в виде соответственно билинейных и линейной форм 
узловых параметров: 


мот вот Г от 
а(и у) =у Ах, Би м)=Уу Вх, с(иу)=Уу Сх, Шу) =ЕУЬ, (6) 
где х и у-— векторы узловых параметров соответственно функций и" и у', а элементы 


матриц жесткости (А), масс (В), демпфирования (С) и вектора нагрузок (Ъ) вычисля- 
ются по формулам (т, / = 1,2, ..., и) 


ау а(фь, фу, Ь, = Ь(Фь фу), су= с(фь фу), В: = Ку, ф). (7) 

Учитывая, что базисные функции подпространства МКЭ Ц’ можно выбрать так, 
чтобы они были отличными от нуля только на нескольких конечных элементах, 
структура матриц из (6) является разреженной и определяется нумерацией узловых 
параметров, которая, в свою очередь, зависит от нумерации узлов конечноэлементной 
сетки (матрица называется разреженной, если число ее элементов, отличных от нуля, 
намного меньше общего их числа — и"). 

Таким образом, учитывая в (6) произвольность вектора у, разрешающие дискрет- 
ные задачи имеют вид [9]: 
— системы линейных алгебраических уравнений (СЛАУ) 


Ах=Ь (8) 
для статической задачи (1); 
— задачи с начальными условиями 


Вх" (0 + Ск + АХО =Ы®, хи) =х®, хи) =х® (9) 


для динамической задачи (2); 
— обобщенной алгебраической проблемы собственных значений 


Ах = ААВх (10) 


для задачи о собственных колебаниях (3). 


Достоверность приближенных решений 


Когда речь идет о достоверности решения, то имеется в виду, что исследуется 
достоверность приближенных решений математических задач (1) -— (3). Полная 
погрешность приближенного решения учитывает влияние погрешностей исходных 
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данных, погрешностей дискретизации и формирования данных соответствующей 
дискретной задачи (8), (9) или (10) и погрешностей компьютерного решения этой 
дискретной задачи. 

Особое внимание следует уделить изучению влияния приближенных исходных 
данных прикладной задачи. Ведь в реальных задачах расчета прочности конструкций 
почти все исходные данные — константы, которыми задаются механические свойства 
материала конструкции (например, модуль Юнга), геометрические размеры, нагрузки и 
тому подобное, как правило, задаются с погрешностями. Для учета влияния погрешнос- 
тей в исходных данных на решения задач теории упругости (поведение наследственной 
погрешности) необходимо исследовать устойчивость решения задачи по коэффициен- 
там уравнений и правой части. Задача называется коэффициентно устойчивой, если при 
малом возмущении коэффициентов задачи ее решение также меняется мало. В [10] 
рассмотрен вопрос коэффициентной устойчивости для первой основной задачи 
теории упругости (когда на границе области заданы напряжения). Описанный в этой 
работе метод исследования коэффициентной устойчивости первой основной краевой 
задачи теории упругости можно использовать для других задач теории упругости и 
для исследования устойчивости по правой части. 

Для МКЭ в перемещениях известны условия сходимости и оценки погрешнос- 
ти решения [3], [8], [9]. Необходимыми условиями сходимости решения, полученного 
с помощью МКЭ, являются выполнение условий существования и единственности 
решения, линейная независимость и полнота системы базисных функций, а также их 
совместность (конформность) либо выполнение условий, компенсирующих несов- 
местность. Линейная независимость следует из равенств (4). Совместность означает, 
что все базисные функции являются возможными перемещениями, т.е. принадлежат 
пространству (0%. 

Известны легко проверяемые условия, позволяющие установить полноту базисных 
функций, их совместность или выполнение условий, компенсирующих несовмест- 
ность. Эти условия имеют вид равенств, которым должны удовлетворять базисные 
функции на каждом конечном элементе. 

При дополнительных условиях на гладкость решений задач (1) — (3) точность 
метода конечных элементов характеризуется априорной оценкой |и-и"|| < С”, где 
й — наибольший из линейных размеров конечных элементов, А — величина, зависящая 
от степени кусочно-полиномиальных базисных функций, которую называют поряд- 
ком сходимости, С\ — постоянная, зависящая от гладкости решения задачи. 

К сожалению, приведенная оценка имеет асимптотический характер, выполняясь 
для достаточно малого й, и не является конструктивной из-за неопределенной кон- 
станты С1. Однако эта оценка позволяет выбрать для решения конкретной задачи 
конечноэлементный базис. Этот выбор обусловливается математической моделью, 
гладкостью решения исходной задачи, необходимой точностью решения, скоростью 
сходимости МКЭ, а также свойствами разрешающих задач (структурой и порядком 
матриц, числом обусловленности матрицы), точностью задания исходных данных. 
Следует также учитывать, что сколь угодно высокая (допускаемая свойствами задачи) 
точность конечноэлементного решения, как правило, является дорогой, да и не всегда 
необходимой. 

Таким образом, теоретические основы МКЭ позволяют не только исследовать 
корректность применения известных конечных элементов, но также разработать прин- 
ципы конструирования новых совместных или несовместных элементов и получить 
для них оценки погрешности [11]. 
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Существенное влияние на точность решения дискретных задач (8) — (10) могут 
иметь погрешности, которые возникают в результате использования квадратурных 
формул при вычислении (7) элементов матриц (жесткостей, масс, демпфирования) и 
вектора нагрузок соответствующей задачи. Чаще всего используются квадратурные 
формулы Гаусса, которые обеспечивают наивысшую алгебраическую точность при 
одном и том же количестве узлов квадратур и для которых справедлива следующая 
оценка остаточного члена [12]: 


- (п!) а}" 
ри 9 


где [2,6] — отрезок интегрирования, п — количество узлов квадратуры. Важно так 
выбирать количество узлов для численного интегрирования, чтобы его точность по й 
была не меньшей точности МКЭ. Кроме этого нужно учитывать, что постоянная С> 
равна максимальному значению на отрезке [а,Б] абсолютной величины производной 
(2п)-го порядка подынтегральной функции. Поэтому абсолютное значение погрешнос- 
ти интегрирования может в некоторых случаях существенно превысить абсолютное 
значение погрешности МКЭ, например, если коэффициенты исходных задач значитель- 
но изменяются в пределах какого-то конечного элемента. В таком случае следует 
увеличивать количество узлов интегрирования (на | — 2 в каждом координатном 
направлении). 

Погрешность компьютерного решения Хх относительно точного решения х 
соответствующей дискретной задачи (8), (9) или (10) зависит от точности задания 
исходных данных, метода и алгоритма решения этой дискретной задачи, а также 
погрешностей округления. Так, полная погрешность решения СЛАУ (8) в случае поло- 
жительно определенной матрицы жесткости А оценивается следующим образом [13] 


хх] < (4) ЕЙ 
У о т 


ый 


где г=Ь- АХ - невязка, й(А) — число обусловленности матрицы А или его оценка, 
&4 = &8 = тазйеру — относительные погрешности задания в компьютере элементов 
матрицы и правых частей системы, тазйеру — наименьшее положительное число с 
плавающей запятой в компьютере, для которого выполняется соотношение 
1 + тазйерх > 1. 

В работе [14] рассмотрено получение с помощью МКЭ достоверного прибли- 
жения к единственному решению задачи на подпространстве (например, в случае 
задания на границе области напряжений). 

После получения решения задачи необходимо провести апостериорный анализ 
результатов с точки зрения удовлетворения тем гипотезам, которые были использо- 
ваны при построении математической модели задачи, а также требованиям точности 
решения. Если решение удовлетворяет всем поставленным требованиям, то можно 
считать, что этому решению можно доверять и использовать для проектирования и 
тому подобное. 


Формирование и решение дискретных задач 


С точки зрения компьютерного решения разрешающих задач (8) — (10) сущест- 
венное значение имеет ряд их особенностей: 
— порядок матриц разрешающих дискретных задач составляет О(10^) — О(10”); 
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— матрицы дискретных задач являются симметричными, положительно определен- 
ными или положительно полуопределенными; 

— матрицы дискретных задач имеют разреженную структуру (ленточную, профиль- 
ную, блочную и т.п.). 

С целью минимизации количества арифметических операций при решении систем 
линейных алгебраических уравнений, а также требуемой оперативной и внешней 
памяти, учитывая разреженную структуру матрицы жесткости А задач (8) - (10), в 
большинстве случаев проводится переупорядочение неизвестных задачи. В зависи- 
мости от критерия минимизации — ширины ленты матрицы, ее профиля, общего 
количества арифметических операций при треугольном разложении матрицы и т. д. — 
используются различные методы переупорядочения [14], например, обратный алгоритм 
Катхилла-Макки, фактор деревьев, минимальной степени. Конкретных рекоменда- 
ций для выбора метода упорядочения дать нельзя, так как эффективность того или 
иного алгоритма существенно зависит от первоначальной структуры конкретной 
матрицы А. 

Формирование дискретных задач. Важным преимуществом МКЭ является то, 
что элементы матриц и векторов правых частей задач (8) — (10) получаются сумми- 
рованием соответствующих элементов матриц (жесткости, масс, демпфирования) и 
векторов нагрузок, построенных для отдельных конечных элементов. Такое свойство 
матриц и векторов дискретных задач МКЭ позволяет эффективно распараллелить 
процесс формирования этих матриц и векторов. 

Например, (метод декомпозиции области) область может быть разбита на при- 
близительно равные подобласти (по числу используемых процессов на основе требо- 
ваний их равномерной загрузки) и для конечных элементов, содержащихся в этих 
подобластях, параллельно и независимо друг от друга формируются соответствую- 
щие матрицы и векторы. После этого с использованием обменов между процессами 
формируются глобальные матрицы и векторы соответствующей дискретной задачи 
(8) — (10). При этом глобальные матрицы и векторы могут быть сформированы или в 
файле (файлах) на внешнем носителе, или распределенными между процессами в 
соответствии с требованиями метода, который будет использоваться для решения 
дискретной задачи. 

Если нумерация узлов и узловых параметров выполняется так, что узлы, распо- 
ложенные на общих границах подобластей (внутренние граничные узлы), получают 
последние номера при сплошной нумерации узлов всей сеточной области, то такой 
подход позволяет сформировать блочно-диагональную глобальную матрицу с окайм- 
лением. В этом случае эффективность распараллеливания близка к единице, так как 
обмены между процессами необходимы только при обработке конечных элементов, 
которые содержат внутренние граничные узлы. 

Во многих параллельных алгоритмах решения задач (8) — (10) используется 
циклическое распределение строк или столбцов матрицы между параллельными про- 
цессами. Тогда циклический алгоритм формирования глобальных матриц МКЭ опи- 
сывается следующим образом: каждый процесс выполняет в цикле для всех групп 
конечных элементов такие действия: 

(1) вычисление всех локальных матриц и векторов подгруппы конечных элемен- 
тов, которая обрабатывается данным процессом на данном шаге цикла; 

(11) в цикле для 1 = 0, 1..., р-1 проводится (р — количество процессов): 
— рассылка из процесса с логическим номером 1 всех вычисленных локальных мат- 
риц и векторов; 
— ассемблирование элементов глобальных матриц и векторов, которые распреде- 
ляются данному процессу, используя локальные матрицы и векторы текущей группы 
конечных элементов, которые сформированы процессом с логическим номером /. 
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Этот алгоритм является наиболее эффективным при циклических схемах рас- 
пределения матриц и векторов между процессами. 

Решение систем линейных алгебраических уравнений. Для решения СЛАУ (8) — 
исходной или переупорядоченной — на МПМО-компьютере могут быть использованы 
параллельные алгоритмы для исследования и решения СЛАУ с симметричными 
матрицами, использующие [ОЁ/-разложение метода Холецкого [15]. В зависимости 
от ширины и заполненности ленты матрицы используются: 

— блочные параллельные алгоритмы для узких ленточных матриц и блочно-диаго- 
нальных матриц с окаймлением, 

— одномерные блочные циклические параллельные алгоритмы для исследования и 
решения СЛАУ с ленточными или с матрицами профильными матрицами. 

При выборе параллельного алгоритма необходимо учитывать его эффективность. 
Например, коэффициент эффективности блочного параллельного алгоритма ГОГ/-разло- 
жение метода Холецкого для узкой ленточной симметричной матрицы оценивается так: 


1 
о инь) 
} 8т 16п 3 2 


а одномерного блочного циклического параллельного алгоритма ГОГ/-разложение 
метода Холецкого для ленточной симметричной матрицы так: 


_ | г 
Е, [1+2 0+0 р ара, 
т т 
2 


где т — полуширина ленты матрицы, 5 — количество строк в блоке, т, =, +—т., 
т 


1 
т. ЕТ, +—®,, 1,1 — отношения времен, необходимых соответственно для установ- 
57 
ления связей между двумя процессами и для обмена одним машинным словом между 
двумя процессами, к среднему времени выполнения одной арифметической операции. 
Решение задачи с начальными условиями. В полудискретном МКЭ прибли- 
женное решение ищется в виде (5), где коэффициенты х; являются функциями 
времени Е. В результате получается система обыкновенных дифференциальных 
уравнений второго порядка с начальными условиями (9), где: х(0), хх векторы 
0 0 | | 
с элементами хх0, х( )= Ё: [® )} х (0 = й: ( )}. 
Система (9) решается методом разложения по формам собственных колебаний. 


Если А, 2. [№ ТВ" =, &=Ь2....,п | — решение алгебраической проблемы собственных 
значений (10), то, полагая в (9) х(г)= уз Ух (()2! ‚ получим вместо (10) (при В-ортого- 
= 


нальности векторов 2 и определенных предположениях относительно матрицы демпфи- 
рования С) систему, которая распадается на независимые уравнения относительно УК®: 


У: (1) +264, у, (1) + у, (В =Р,(®, 1>4, 


(0) (0 


(11) 
Ук (1) = Ур, Ук) = У, 


Е т ИЕ т 
где м =А;*, 0<&<1, Р@Феаь БО, у =х® Ва, уд хо ВЕ, 2, дп. 
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Решения задач (11) имеют вид: 


(1) (0) 
щи Ув № У, 0) 


У, 


Е + УФ со5, Е |+ 


у, (#) =е 


= -| Р,(т) 


й 
9), ^° 


—в ких ((-Т) 


тб), (1— тат, 


где ®, =®,./1-&, . Анализ этих решений свидетельствует о том, что существенный 


вклад в решение х(1) задачи (9) вносят лишь около 10 составляющих, соответствующих 
минимальным собственным значениям. Поэтому соответствующая задача на собственные 
значения вида (10) может решаться методом итераций на подпространстве [16]. 

Существует ряд случаев нагружений, когда возможно точное вычисление решения 
задачи (11). Например, для ветровой, сейсмической или гармонической нагрузок. 
В остальных случаях решение у) находится численно, например, методом конечных 
разностей по схеме Ньюмарка. Также задачи Коши (11) на МИМО-компьютере могут 
решается методом Рунге-Кутта четвертого порядка, ориентированным на решение 
систем обыкновенных дифференциальных уравнений второго порядка. 

Решение обобщенной алгебраической проблемы собственных значений. Решение 
на ММО-компьютере обобщенной алгебраической проблемы собственных значений (10) 
с ленточными или профильными симметричными матрицами (исходными или переупоря- 
доченными) может быть получено с помощью параллельных алгоритмов метода итераций 
на подпространстве, которые, как и в случае решения СЛАУ (8), используют одномерное 
блочно-циклическое распределение данных между процессами [17]. 


Решение задач расчета строительных 
конструкций на Инпарком 


Как отмечалось выше, для исследования и решения задач (1) — (3) может быть ис- 
пользован программный комплекс (ПК) по решению некоторых задач расчета прочности 
строительных конструкций на многопроцессорных рабочих станциях Инпарком. Ниже 
представлена технологическая схема этого программного комплекса (табл. 1). 


Таблица 1 — Технологическая схема программного комплекса Инпарком 


Обработка 
Формирование 
Создание анны Решение результатов Анализ 
КЭ-модели ОЕ разрешающей решения результатов 
конструкции р и задачи разрешающей расчета 
задачи 
интерактивно 
на языке пред- автоматически автоматически автоматически интерактивно 
метной области 
препроцессор препроцессор пра Ъ препроцессор препроцессор 
ПК ЛИРА ПК ЛИРА 9.2 р ПК ЛИРА 9.2 ПК ЛИРА 
локально локально 
на хост- на хост- 
локально локально локально 
компьютере компьютере 
на хост- на на хост- 
Инпарком или Инпарком или 
компьютере обрабатывающем компьютере 
удаленно, И удаленно, 
нпарком блоке Инпарком Инпарком 
используя используя 
Интернет Интернет 
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В этом случае сформированные на графической станции препроцессором про- 
граммного комплекса исходные данные соответствующей дискретной задачи (8) — (10) 
записываются в бинарные файлы, используя принятые в ПК ЛИРА форматы данных. 
После этого с помощью специального интерфейса на обрабатывающие узлы компьютера 
Инпарком загружается параллельная программа решения соответствующей дискретной 
задачи, которая: 

— анализируя параметры задачи, определяет алгоритм решения задачи и параметры, 
необходимые для решения; 

— считывает исходные данные из файлов и распределяет их между процессами в 
соответствии с выбранным параллельным алгоритмом решения задачи; 

— исследует и решает задачу, обращаясь к соответствующим функциональным моду- 
лям библиотеки шраг!Ь [4], [5]; 

— сохраняет результаты исследования и решения задачи в файлах для дальнейшей 
обработки и использования постпроцессором программного комплекса. 

Результатами статического расчета напряженно-деформированного состояния 
являются перемещения узлов схемы и усилия (напряжения) в сечениях элементов. 
Результатами динамического расчета — решения АПСЗ (10) — являются периоды, час- 
тоты и формы собственных колебаний для каждого тона. Результатами решения 
задачи Коши (11) являются значения перемещений ух; = ук[) в узлах. 

Программным комплексом по решению некоторых задач расчета прочности 
строительных конструкций на рабочих станциях Инпарком проводилось решение ряда 
тестовых статических задач, для которых разрешающими являются СЛАУ (8). 
В табл. 2 приведены времена решения СЛАУ тестовых задач на рабочих станциях 
Инпарком, а также ориентировочные времена решения этих задач на графической 
станции компьютера Инпарком программным комплексом ЛИРА 9.2. 

Следует обратить внимание, что проводимые функциональными модулями пра 
исследования увеличивают общее время решения задачи, как правило, не более, чем на 
10 %. Однако такие исследования (особенно вычисление оценки числа обусловленности 
матрицы СЛАУ и сопоставление этого значения с погрешностями исходных данных 
задачи при вычислении оценки наследственной погрешности) позволяют оценить досто- 
верность полученного решения и, что особенно важно, внести коррективы в постановку 
задачи при невозможности гарантировать достоверность решения задачи. 


Таблица 2 — Времена решения СЛАУ тестовых задач 


№ Порядок и - ое. Ом ЕеНИЯ ИЕ) Ононение 

п/п системы (макс.) | ность а Инпарком ПК Лира а 
вЫ процессов т решения 
64 0,200 9,55 

| 444 44 21° 1 
29 :о й 96 0,183 я 10,42 
16 0,117 6,51 

2 44 436 37 580 2% - : 
- 80 0,083 Бы 9,12 
3 189 956 22 585 9% 16 1,937 35,58 18,37 
4 283 031 19 530 7% 16 1,222 16,40 23,42 
16 10,083 5,38 

| 4242 о 4.24 
. а ] 2. 96 5,200 сы 10,43 
6 1000 332 1004 85% 16 3,001 41,00 13,66 
7 1200 000 1265 | 100% 16 4,776 102,41 21,36 


Пример расчета на прочность конструкции здания (рис. 1) с помощью програм- 
много комплекса по решению некоторых задач расчета прочности строительных 
конструкций на рабочих станциях Инпарком. 
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Рисунок 1 — Конструкция здания Рисунок 2 — Структура 
для примера расчета на прочность переупорядоченной матрицы 


Конструкция разбита на 972 808 конечных элементов. Образованная в резуль- 
тате разбиения конечно-элементная сетка содержит 895 302 узла. 

В результате дискретизации получена разрешающая задача (8) — система линейных 
алгебраических уравнений, порядок которой равен 5 371 727. После оптимизации струк- 
туры матрицы (переупорядочения) методом минимальной степени плотность (заполнен- 
ность ленты матрицы ненулевыми элементами) переупорядоченной структуры составила 
менее 1 %, хотя полуширина ленты стала равной 3 174 975. Структура переупорядочен- 
ной матрицы представлена на рис. 2. Размер файла для хранения информации о структуре 
матрицы и значений ее ненулевых элементов составляет около 76 ОБ. 

Задача была сформирована и решалась с помощью программного комплекса на 
рабочих станциях Инпарком. Для решения разрешающей системы (8) использовались 
функциональные модули библиотеки шра!Ъ. Расчеты проводились для одного случая 
нагружения. Для расчетов использовалось различное количество процессоров (и про- 
цессов) рабочих станций Инпарком-64 и Инпарком 128. В табл. 3 приведены времена 
решения дискретной задачи. 


Таблица 3 — Времена решения дискретной задачи 


Компьютер Инпарком-128 Инпарком-64 а 
с 32 24 20 16 12 | 
процессоров 
Количество 

. процессов 128 96 80 64 48 > 1 
ремя решения . : . р | часов 

СЛАУ 6:09 7:27 8:18 9:24 10:22 (оценка) 
Выводы 


Проведенные авторами исследования и апробация программного комплекса по 
решению некоторых задач расчета прочности строительных конструкций на рабочих 
станциях Инпарком показали перспективность создания программных средств по реше- 
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нию прикладных задач на параллельных компьютерах путем встраивания параллельных 
программ по решению отдельных подзадач в существующие программные комплексы 
для решения прикладных задач в различных отраслях науки и техники. В качестве таких 
параллельных программ могут использоваться многие программы библиотеки шра® 
или ее аналогов на других компьютерах с параллельной обработкой информации. 

Также является перспективным расширение круга решаемых задач с помощью про- 
граммного комплекса по решению некоторых задач расчета прочности строительных 
конструкций на рабочих станциях Инпарком и распараллеливание формирования данных 
дискретной задачи и обработки результатов решения этой задачи, используя предложен- 
ные в настоящей работе алгоритмы. 
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О.М. Хи, В.В. Полянко, О.В. Попов, О.В. Рудич 

Розв’язування задач розрахунку мщност! конструкций на МИМО-комп’ютер! 

Розглядаються блочн! алгоритми дослдження та розв’язування задач лийно! алгебри з розрдженими 
(вузькими страчковими, з обрамленням 1т.п.) симетричними матрицями на комп’ютерах МПМО-архтектури. 
Дослджуеться ефективнсть даних алгоритм. Наводяться деяк! результати чисельних експериментив на 
М1МО-комп’ютери. 


А.М. Юитсй, Г.Г. Рорзапко, А.Г. Роротъ, О.Т. Ки@сй 

'Тве Зо!ушо о? РгоШепа$ оЁ Эгисига! Апа!у$15 оп МИМО-сотрщег 

Те зо[уше оЁ ргоетл$ оЁ зтас@га] апа]уз1$ оп рагаПе| согаршег$ аге 4еай у, шста@те а]еогИиа$ оЁ 
Гог, шуезисайпе ап зо]уше оф зузепа$ оЁ Фе зе пе едиайоп$ гесегуеа Фигте @1зсгейтаноп Бу Фе 
ВНоНе е@]етеп тео. Зоте гези 5 оЁ питега| 1е515 сагле4 ой оп МПМО-сотршег аге э1уеп. 
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